Magnetic domain wall motion in a nanowire: depinning and creep 
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The domain wall motion in a magnetic nanowire is examined theoretically in the regime where 
the domain wall driving force is weak and its competition against disorders is assisted by thermal 
agitations. Two types of driving forces are considered; magnetic field and current. While the field 
induces the domain wall motion through the Zeeman energy, the current induces the domain wall 
motion by generating the spin transfer torque, of which effects in this regime remain controversial. 
The spin transfer torque has two mutually orthogonal vector components, the adiabatic spin transfer 
torque and the nonadiabatic spin transfer torque. We investigate separate effects of the two compo- 
nents on the domain wall depinning rate in one-dimensional systems and on the domain wall creep 
velocity in two-dimensional systems, both below the Walker breakdown threshold. In addition to 
the leading order contribution coming from the field and/or the nonadiabatic spin transfer torque, 
we find that the adiabatic spin transfer torque generates corrections, which can be of relevance for 
an unambiguous analysis of experimental results. For instance, it is demonstrated that the neglect 
of the corrections in experimental analysis may lead to incorrect evaluation of the nonadiabaticity 
parameter. Effects of the Rashba spin-orbit coupling on the domain wall motion are also analyzed. 



I. INTRODUCTION 

A magnetic domain wall (DW) in a ferromagnetic 
nanowire is an important subject in spintronics. A new 
type of logic device is proposecU- based on the DW dy- 
namics and a DW-based memory is also proposed^, which 
may have merits such as nonvolatility, high speed, high 
density, and low power consumption. 

The dynamics of a DW varies considerably depending 
on the relative strength of DW driving forces (such as a 
magnetic field and a current) with respect to disorders, 
which tend to suppress the DW motion. If the forces are 
sufficiently strong or the disorders are sufficiently wcak^, 
the DW dynamics does not deviate much from the ideal 
dynamics in the absence of disorders. While some exper- 
iments^^— are estimated to be in this regime, many other 
experiments^^— appear to be in the regime where the 
disorders are important. It is thus desired to understand 
the DW dynamics in the weak driving force regime where 
the competition between the DW driving forces and the 
disorders is significant. 

The DW motion in the weak driving force regime is 
an important example in the field of driven interfaces. 
The study on driven interfaces has a long historyi^ and 
addresses many physical systems such as surface growth 
of a crystali^, vortex line motion in high temperature 
superconductorsi^, and fluid propagation in porous me- 
diapi^. Through a long series of theoretical worksi^ii^iii, 
a simple picture has emerged; the interface motion be- 
comes collective and the collective length scale Lcor^, 
which characterizes the length scale of collectively mov- 
ing interface segments, diverges in the weak driving force 
limit. Due to the divergence of icoi, the interface has 
to overcome an increasingly larger energy barrier as the 
driving force becomes weaker, with the energy barrier Eb 
as a function of the driving force / diverging as a power 



law, Eb oc (/i > 0). Interestingly the creep exponent 
/i is universal in the sense that its value does not change 
continuously with variations of system details and is af- 
fected only by a small number of key features such as the 
system dimensionality. Systems with the same exponent 
are said to be in the same universality class. 

This prediction has been unambiguously confirmed 
for the field-driven DW motion in metallic ferromag- 
nets^, where the DW velocity v is proportional to 
exp(— K_ff~^/fcBr). Here ks is the Boltzmann con- 
stant, T is the temperature, H is the magnetic field 
strength, and k is a constant. Note that this behav- 
ior of w is a combined result of the power law scaling of 
the energy barrier Eb = kH~^ and the Arrhenius lawi^ 
V (X exp(— E^b/ZcbT). The creep exponent ^ is found to 
be sa 0.25, which agrees with the theoretically predicted 
value 1/4 in two-dimensional (2D) systems^. 

A pioneering experiment)^ revealed interesting twists. 
For nanowires made of a ferromagnetic semiconductor 
(Ga,Mn)As, the energy barrier for the field-driven DW 
motion was found to scale as H^^, where fi w 1.2 in- 
stead of 1/4. This difference was attributed to the differ- 
ent nature of disorders; while disorder potential energy 
is short-range correlated in metallic fcrromagnets, it was 
argued that in ferromagnetic semiconductors, disorder 
force is short-range correlated. Since the disorder poten- 
tial energy is obtained by integrating the disorder force, 
it implies that the disorder potential energy is then long- 
range correlatedii. For such cases, it is knownii that 
the nature of the correlation along the DW segments is 
modified and the value of fi indeed changes. 

Another interesting twist of the experiment^ is that 
for the current-driven DW motion, the effective energy 
barrier was reported to scale as where J is the 

current density and fi w 0.33 rather than 1/4 or 1.2. 
Thus two different creep exponent values (1.2 and 0.33) 



2 



were obtained from the same material, implying that the 
current-driven DW motion is qualitatively different from 
the field-driven DW motion. 

It is believed that the current induces the DW mo- 
tion in a nanowire by generating the spin transfer torque 
(STT). The STT has two mutually orthogonal vector 
components, the adiabatic STT and the nonadiabatic 
STT— ii^. The nonadibatic STT has similar properties as 
the magnetic field while the adiabatic STT has very dif- 
ferent properties. Thus the experimental result^ implies 
that the nonadiabatic STT cannot be the main driving 
force of the DW motion. In fact the exponent fi w 0.33 
has been interpreted^ as an indication that the current- 
driven DW motion is mainly due to the adiabatic STT. 

This interpretation is at odds, though not contradic- 
tory, with other results. In metallic ferromagnets, the 
onset of the adiabatic-STT-driven DW motion is esti- 
mated2£i2i to occur at the current density of ^ 10^ 
A/cm^, which is unendurably high for most experimental 
systems. Thus the DW motion realized at lower current 
densities arc usually attributed to the nonadibatic STT. 

This situation strongly motivates experimental^^ and 
theoretical^^"— studies of the current-driven DW motion 
in metallic ferromagnets. This paper aims at theoretical 
explorations of this issue based on the observation that 
the DW anisotropy, characterizing the energy cost asso- 
ciated with the change in the tilting angle of the mag- 
netization inside a DW , is orders of magnitude larger 
in metallic ferromagnets than in ferromagnetic semicon- 
ductors. Since the DW anisotropy tends to suppress 
variations of the tilting angle, we assume that the DW 
creep motion in metallic ferromagnets exhibits the below- 
the-Walker-brcakdown-like behavior in the sense that the 
amplitude of the tilting angle variations during the creep 
motion stays much smaller than 2tt. For the field-driven 
DW creep motion, this assumption is experimentally sup- 
ported since the experimental value ^ 0.25 of the creep 
exponent agrees with the prediction 1/4 of the thcoryii, 
in which the tilting angle dynamics is completely sup- 
pressed. For the current-driven DW creep motion, the 
assumption requires an experimental confirmation. A re- 
cent experiments^ reports the purely current-driven DW 
creep motion in metallic ferromagnets. For ferromagnetic 
semiconductors, in contrast, it appears that the assump- 
tion may not be valid. For the current-driven DW creep 
motion, it was argued^ that each thermally-assisted tun- 
neling event overcomes the energy barrier generated by 
the DW anisotropy, implying that each tunneling event 
is accompanied by the tilting angle change by ^ tt. 

The paper is structured as follows. In Sec. |lll we dis- 
cuss the DW depinning from a single potential well in 
one-dimensional (ID) systems. Analysis of this relatively 
simple problem clearly illustrates separate roles of the 
magnetic field, the adiabatic STT, and the nonadiabatic 
STT on the thermally-assisted tunneling of a DW. It also 
allows one to identify relevant factors affecting the tun- 
neling, which therefore should be included in the analysis. 
In this sense, Sec.|ll]is pedagogical. Nevertheless predic- 



tions in Sec. |TT] can be tested in real experiments since 
a DW exhibits the ID dynamics^i when L^oi becomes 
larger than both the thickness and width of a nanowire. 
In particular, it is predicted that when the depinning 
rate is used as a tool to evaluate the nonadiabaticity pa- 
ramete r^^'^^ , characterizing the strength of the nonadi- 
batic STT, it may lead to incorrect values if disorders in 
a nanowire have certain features. In Sec. IIIIl the DW 
creep motion in 2D systems is analyzed. Separate roles 
of the magnetic field, the adiabatic STT, and the nonadi- 
abatic STT on the creep motion are clarified. In addition 
to the leading order contribution to the creep motion in 
the vanishing DW driving force limit, next leading or- 
der contributions are also obtained. Although the next 
leading order contributions are irrelevant as far as the 
theoretical determination of the creep exponent and the 
universality class is concerned, they may nevertheless be 
relevant in experimental determination of the creep expo- 
nent since experiments are always performed at small but 
finite driving force strength. At the end of both Sees. HIl 
and lllli effects of the Rashba spin-orbit coupling (RSOC) 
arc discussed. The emergence of the RSOC in ferromag- 
netic nanowircs is recently demonstrated^^. Section IIVI 
concludes this paper. 



II. DW DEPINNING IN ID 

When both the thickness and the width of a mag- 
netic nanowire are sufficiently smaller than the collective 
length Lcoi, the system reduces to a ID problem and the 
configuration of a DW can be described by two variables, 
the DW position q and the tilting angle -0. This Section 
examines the DW depinning from a potential well in this 
ID regime. 



A. Effective energy 

In the ID regime, the response of the DW collective co- 
ordinates {q,tp) to an external magnetic field H and/or 
an electric current of density J is described by the fol- 
lowing equations, 
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2MsnX d^jj ' 



(1) 
(2) 



where a is the Gilbert damping parameter, A is the DW 
width, 7o is the gyromagnetic ratio, Ms is the satura- 
tion magnetization, D, is the cross-sectional area of a 
nanowire, V{q, ip) is the DW potential energy, and the 
dimensionless parameter /3 is the nonadiabaticity coeffi- 
cient^ii^ representing the strength of the nonadiabatic 
STT. X = hP/2XeMs (< 0) is a constant with the di- 
mension H/ J, P is the spin-polarization of the current, 
and h is the Planck constant. In Eqs. ([T]) and the 
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Fig 1: (Color online) Schematic plot of the energy landscape 
E{q,tp) as a function of q and '0. (a) The activation path 
(black solid arrow) of the thermally-assisted transition for 
H ^ and J = is shown schematically from one local mini- 
mum A [(qajiPa) ~ {qcOjipGo)] to another local minimum C 
{qc,ipc) through the saddle point B [(i7s,'(/'b) = iqso,tpso)]- 
For simplicity, ^pco = tpso ~ V'c = is assumed, and 
E{qGo,il^Go) > E{qc,tpc) is also assumed in this plot, (b) 
When J is turned on, the activation path (black solid arrow) 
is deformed to A'-B'-C. Here the coordinates of A' and B' 
correspond to (qcil^a) and {qs,i^s) in the text, respectively. 
Note that the path is now curved due to the adiabatic STT 
when Uq'^ ^ ^s'^- 



sign convention of J is chosen in such a way that positive 
J drives the DW towards the positive q direction. On 
the other hand, the sign convention of H should depend 
on the types of the DW; In nanowires with the in-plane 
anisotropy, for instance, opposite signs should be adopted 
for the head-to-head and tail-to-tail DWs. Below for sim- 
plicity, we consider one particular sign only, and assume 
that the positive H tends to drive the DW towards the 
positive q direction. 

The DW potential energy V{q^ ip) consists of the DW 
anisotropy energy 2{lXKd sin^ ip and a disorder poten- 
tial energy, where Kd represents the strength of the DW 
anisotropy. Here tp is defined in a way that tp — for 
the tilting angle preferred by the DW anisotropy. When 
the disorder potential energy depends only on q, Eqs. (H]) 
and become equivalent to Eqs. (3) and (4) in Ref. |23. 
In general, however, the disorder potential energy may 
also depend on -0. 

Such ip dependence of V may arise in various ways. 
For instance, the value of Ms may fluctuate from posi- 
tion to position. RecaUing that Kd depends^i on Ms, 
Kd can then be decomposed into its spatial average part 
and the fluctuating part 5Kd{q). The fluctuating part 
of the DW anisotropy energy [oc 5Kd{q) s\v? ip] may be 
absorbed to V to generate its ip dependence. Similar 
dependence may arise from the position-to-position fluc- 
tuation of f2. In these types of disorder, the preferred 
tilting angle remains unaffected and only the strength of 
the DW anisotropy fluctuates. Some types of disorder 
may generate opposite effects. In a magnetic nanowire 
with the perpendicular magnetic anisotropy (PMA), the 
interface between the magnetic layer and the neighboring 
layer plays important roles for the anisotropy. When the 
interface is not perfectly flat and becomes rough^^Ti^, the 
preferred anisotropy direction fluctuates from position to 



position. In this case, the preferred tilting angle fluctu- 
ates while the strength of the DW anisotropy may not 
fluctuate. 

Below we consider this general situation, in which the 
disorder potential energy depends both on q and ip. In 
Ref. the ip dependence of the disorder potential en- 
ergy is included in its initial formulation but ignored 
when the depinning rate is calculated. We demonstrate 
below that the ip dependence of V generates interesting 
consequences. 

Based on the Lagrangian formulation, Eqs. ([1]) and ^ 
may be considered as the Lagrange's equations of the 
Lagrangian L and the dissipation function F, 

L=^{qiP-qi;)-V{q,iP) (3) 
7o 

+ 2Msnq{H - PxJ) - 2MsnXijxJ, 

F^a^ie + X^'). (4) 
70 A 

The Lagrangian in Eq. (jS]) is then transformed to the 
Hamiltonian i.e. the effective energy function E, 

Eiq,lP)^V{q,lP)-2MsnqiH -pxJ) (5) 
+ 2MsVLXipxJ- 

Here we have used the term effective energy since E is 
not a single valued function^ in the sense that E{q, ip) ^ 
E{q,ip + 2TT) although {q,ip) and {q,ip + 2'K) represent the 
same magnetic configuration. Thus some care should be 
exercised when Eq. ([5]) is used to analyze the DW dy- 
namics above the Walker breakdown threshold, where ip 
changes more than 2tt. Below the Walker breakdown 
threshold, on the other hand, the dynamics of ip is con- 
fined to a value range narrower than 27r and E{q, ip) can 
be regarded as a single valued function. 

B. Effective energy barrier 

Figure [T] shows schematically the energy profile, to 
which the DW is subject. The DW has to overcome 
an energy barrier to get depinned from a given poten- 
tial well. When the DW driving force {H or J) is small, 
the height of the energy barrier is sufficiently higher than 
the DW energy measured from the bottom of the poten- 
tial well and the DW overcomes the large energy barrier 
by exploiting the thermal agitation. Thus the depinning 
time from potential wells is governed (within the expo- 
nential accuracy) by the energy barrier via the Arrhenius 
law. When the depinning time is much longer than the re- 
laxation time inside potential wells, the energy barrier is 
defined as the difference between the saddle point energy 
and the local ground state energy. One remark is in or- 
der. While the Arrhenius law is based on the fiuctuation- 
dissipation theoreni^^ri^^, the theorem does not generally 
hold when J is finite and the system is thus in nonequilib- 
rium situations. However it has been demonstrated that 
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for small J= and below the Walker breakdown regime^, 
thermal fluctuations still satisfy the theorem, justifying 
the use of the Arrhenius law in this case. 

The energy barrier Eb depends on H and J, and we 
examine this dependence. For H = J = 0, E{q, -iji) re- 
duces to V{q,il)). Let (gGO,'0Go) and (gso,V'so) denote 
respectively the local ground state and saddle point con- 
figurations of '(/')■ Note that wc introduce separate 
parameters ^go and "ipso- Although ■050 — "^go will be 
much smaller than 27r in the regime below the Walker 
breakdown, the difference is nonzero in general due to 
the ^ dependence of the disorder potential energy. To 
examine effects of small H and J, V{q, ip) may be Taylor 
expanded near these configurations; 



(6) 



for (5,-0) near (gccV'Go), and 

V ^ Vo - L^liq ^ qsof + - tpsof, (7) 

for near (gso,'0so)- Here uc/s and vq/s are the 

potential stiffness, and the potential depth Vq amounts 
to the energy barrier height ior H = J = 0. Note that 
uJq in Eq. © and ujg in Eq. (O appear with the opposite 
signs due to the difference between the local ground state 
and saddle point (Fig. [T|). Note also that wc distinguish 
i^G and in order to take account of the tjj dependence 
of the disorder potential energy. 

The driving forces H and J modify the local ground 
state and saddle point configurations to, say, {qcipc) 
and {qs,ips)- For small H and J, the modified configu- 
rations can be determined from 5E = with the aid of 
Eqs. (|6]) and ([7]). One obtains 

qo = 9G0 + Ms^^a" ~ PxJ), (8) 

V-G = ^Go - Ms^Vg^XxJ, (9) 

qs = qso-MsnuJs\H-PxJ), (10) 

= V'so - Ms^vfXxJ, (11) 

The evaluation of the energy barrier Eb = E{qs,ips) ^ 
E^qci^o) is now trivial. One finds. 
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2SqoMsniH - PxJ) 

Mln'uz'\XxJ)\ 



(12) 



where Sqo 

^ "^G 



qso 



qoQ, 



Stpo 



■050 - 0GO, UJ_^_ = 
Lijg~ + ujq~ and I'Z'^ ^ 1^5^ — ^g^- Equation p2p clearly 
shows the effect of H and J on the energy barrier. Among 
the two components of the STT produced by J, the non- 
adiabatic STT (ex (3xJ) in the first and third lines of 
Eq. (fT2)) has the exactly same effect as the magnetic field 
H while the second and the fourth lines of Eq. ((T^ in- 
dicate that the effect of the adiabatic STT (cx Ax J) is 
qualitatively different from the field effect. 



The depinning rate l/r from a potential well is then 
given by 
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Eb{H,J) 



(13) 



where 1 /tq amounts to the attempt frequency and the H 
and J dependence oi Eb are given in Eq. ([T^. 

Recently Kim and Burrowes^^ analyzed the effective 
energy barrier for the purely current-driven DW creep 
motion in ID. Equation (39) in their work indicates that 
J modifies the energy barrier Eb through a linear term 
(oc /3J), and a quadratic term (cx /3^J^), both of which 
arise from the nonadiabatic STT. Our result [first and 
third lines in Eq. (fT2)) ] agrees with this result as far 
as these two terms are concerned. However our result 
predicts that there are another linear term [cx Ax J, the 
second line in Eq. and quadratic term [cx }?x^,P , 

the fourth line in Eq. ((T2|) ]. which arise from the adia- 
batic STT. This difference between our result and Ref.[25l 
stems from the nature of the disorders; In Ref. [2^, the 
calculation oIEb assumed that the disorder contribution 
to Y{iq, ■0) depends only on q and does not depend on 0, 
whereas we consider more realistic situations where the 
disorder contribution depends not only on q but also on 
■0. This dependence on "0 appears in the second and last 
terms in Eq. through the factors (5-0o and vZ^ ■ 

DW depinning experiments^"— are sometimes used 
as a tool to determine the nonadiabaticity parameter /3. 
When the ^0 dependence of the disorder potential energy 
is negligible and thus J^'o = vZ^ ~ 0, one can verify from 
Eq. (fT2|) that Eb depends on H and J through a single 
variable H — PxJ ■ Thus by comparing the "efficiency" 
of H and J in the DW depinning, one can determine 
p. In general, however, the ip dependence of the disorder 
potential energy may not be negligible. In such situations 
and in the limit H,J — > 0, the H and J dependence of 
Eb appears through a different single variable H ~ j3xJ — 
XxJSipo/Sqo thus uncareful experimental evaluation may 
incorrectly identify 



/?' = /5 + A^ 
oqo 



(14) 



as /3. Thus possible -0 dependence of the disorder poten- 
tial energy should be carefully examined for the correct 
evaluation of /3. 

As discussed in Sec. Ill Al the V dependence of V may 
be qualitatively different depending on details of disor- 
ders. When vZ^ ^ but Sipo = 0, the second contribu- 
tion in Eq. (|14[) vanishes, simplifying the experimental 
evaluation of /3. When uZ^ = but ^■00 0> on the 
other hand, the second contribution in Eq. (jl4p may not 
be negligible. A possible way to avoid the incorrect eval- 
uation of /3 in this case is to take an average of /3' for 
multiple potential wells. Since the sign of Sipo is expected 
to fluctuate from potential wells to potential wells, this 
averaging process may be able to remove the second con- 
tribution of /3' proportional to StpQ. By the way, the sign 
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J. From Eq. p2)) . one finds that H* satisfies 
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Fig 2: (Color online) (a) An example with the broken inver- 
sion symmetry. The ferromagnetic layer (FM) is sandwiched 
between two difterent nonmagnetic layers (NMl and NM2), 
so that the inversion symmetry is broken along the y direc- 
tion. When the current is injected along x direction, the 
Rashba spin-orbit coupling makes the magnetization feel as 
if a magnetic field Hk is applied^ along the z direction, (b) 
Schematic plots of the magnetization configurations for the 
Bloch wall (upper plot) and the Neel wall(lower) in nanowires 
with the perpendicular magnetic anisotropy. Solid arrows 
(colored in red) represent local magnetic moments inside of a 
DW. For a Bloch (Neel) wall, i/) represents the angle between 
the magnetization and the positive z(x) axis within the xz 
plane. 



fluctuations of bq^ can be suppressed in this averaging 
process since the depinning to the right ((5go > 0) and to 
the left ((5(7o < 0) are distinguishable in experiments. 

Lastly we compare two contributions [third and fourth 
hues in Eq. (fT2| ]. both of which generate the J-quadratic 
contributions to Eb- They have one important differ- 
ence; the third Une, which arises from the nonadiabatic 
STT, always enhance Eb and thereby lower the depin- 
ning rate while the fourth line, which arises from the 
adiabatic STT, may either increase or decrease Eb since 
vZ^ can be positive or negative depending on the nature 
of disorders. Thus in case that experiments find the J- 
quadratic contribution enhances the depinning rate, it 
impUes that the adiabatic STT makes a larger contri- 
bution to the J-quadratic dependence of Eb than the 
nonadiabatic STT. 



C. Effective magnetic field 

The DW depinning for the purely field-driven case is 
relatively well understood§'Ji. Thus if one can "map" 
general situations with both H and J to the purely field- 
driven case, it may provide a useful way to describe ex- 
perimental results in general situations. The effective 
magnetic field is one way to make this connection. We 
define the effective field H* {H, J) of the DW depinning 
by the relation Eb{H\Q) = Eb{H,J). H*{H,J) can 
be experimentally extracted, for instance, from contour 
plots^ of the DW depinning rate as a function of H and 
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(15) 
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{H - (ixJ) 



(AxJ)-^(AxJ)^ 



Solving Eq. (|15p for H* under the constraint H*{H 
0, J = 0) = leads to 



MsVl MsVl 



1 - 



Sqo - Sq 

. 1/2 



(16) 



^lisqor j ' 

where Sq = qs ^ Qg, Stjj = tjjs — ipG- Here we have 
used the relations, Msfl{H — PxJ) — "^+(^90 ~ Sq) 
and Ms^XxJ = ^-iS^^o — Sip) obtained from Eqs. ([8]), 
dH), (dni), and (HH). Since {Sqo - Sq)/Sqo < 1 and 
{v^[(Sipo)^ — W)^]}/[i^+('5go)^] < 1, one can expand 
the curly bracket in Eq. ([TB| to obtain 



H*{H, J) =H-p'xJ 



2u'iSqo 



(AxJ)' 



(17) 



2uj\Sqo Sqo 



{XXJ){H- P'xJ) 



oiJ'). 



In case the "0 dependence of the disorder potential en- 
ergy is negligible, /?' = /3, i^Z^ = Sipo = 0, and the ef- 
fective field H* reduces to H — PxJ ■ Then the points in 
the (iJ, J) plane with the same depinning rate will form 
straight lines with the slope Px- 

However in more general situations with the V' depen- 
dence of the disorder potential energy, deviations from 
this simple result will occur . When V but Si^a = 0, 
the contour lines of the equi-depinning rate will not be 
straight but instead form parabolas in the (iJ, J) plane 
with the coefficient of the J-quadratic term proportional 
to I'Z^ ■ Note that this quadratic contribution to H* is 
entirely due to the adiabatic STT, while in case of Eb , 
both the adiabatic and nonadiabatic STTs can generate 
the J-quadratic contributions [Eq. (fT2)) ]. In this sense, 
H* allows clearer separation between the adiabatic and 
nonadiabatic STT contributions. When vZ^ = but 
Stjjo ^ 0, the contour lines of the cquiv-depinning rate 
will form straight lines with the modified slope, In 
this case, the value of /?' will fluctuate from potential 
wells to potential wells. 

The above analysis provides experimental procedures 
to determine whether or not the i\) dependence of the dis- 
order potential energy is negligible in a given experiment; 
If the contour lines of the equiv-depinning rate are not 
straight lines, vZZ^ is not zero. If the slope of the lines 
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tangential to the contour lines at the points {H, J = 0) 
fluctuates from potential wells to potential wells, S^q is 
not zero. 



D. Rashba spin-orbit coupling effects 

The special theory of relativity requires the cou- 
pling between the spin and orbital degrees of freedom^. 
Thus the spin-orbit coupling (SOC) is ubiquitous. The 
strength of the SOC however varies considerably from 
systems to systems. It is well known^ that the SOC 
may be considerably enhanced in systems with the bro- 
ken inversion symmetry. The SOC in this case is called 
the Rashba SOC (RSOC). Magnetic systems are not ex- 
ceptions and the RSOC develops in magnetic systems 
with the broken inversion symmetry, as exemplified in a 
recent experiment2£. 

Since the RSOC affects conduction electron spins 
and they in turn interact with the local magnetization 
through the s-d exchange coupling, it also affects the lo- 
cal magnetization. It was reportediS, that a high DW 
velocity can be achieved in magnetic films with the bro- 
ken inversion symmetry. In this subsection, we discuss 
the RSOC effects on the DW depinning. 

When the conduction electron spins are modified by 
the RSOC, according to Ref. H^, the s-d exchange cou- 
pling generates an additional magnetic field acting on the 
local magnetization. Though this is not a real magnetic 
field, it behaves just like a real magnetic field as far as its 
effect on the local magnetization is concerned. When the 
inversion symmetry is broken along y direction and the 
current is injected in x direction[Fig.[21[a)], this magnetic 



field 



^45 



-ffRSOC = — ^^J(£ X y). 



(18) 



where ur is the RSOC constant and is the Bohr 
magneton^. The direction of this field may or may not 
be parallel to the real magnetic field applied to induce the 
DW motion. When it is parallel, its effect is trivial since 
one just needs to replace H hy H + H^goc in all equa- 
tions presented above. When it is not parallel, it may 
induce the current-induced tilting angle jump at strong 
HrsoC: similarly to the chirality switching predicted for 
oblique magnetic field^. For weak Hjisoc, the tilting 
angle jump is unlikely and a separate analysis is required 
to understand its effect on the depinning. 

As a representative example of nonparallel situations, 
we consider a nanowire with the PMA^^i^i^i^i^ along 
the y direction. Then the external magnetic field H 
(along y direction) for the DW motion and Hrsoc (along 
z direction) are mutually orthogonal. In PMA nanowires, 
two types of DWs can exist depending on the width w 
of the nanowiro2i[Fig. [2](b)].; When w is larger than a 
threshold value, a Bloch wall is energetically preferred, 
and when w is smaller than the threshold value, a Neel 
wall is preferred. 



One of primary effects of i?Rsoc is to modify the tj; 
dependence of E[q,'\}j) [Eq. (O], since the energy of the 
system is minimized when the magnetization direction 
at the center of the DW is parallel to i?R,soc- Below 
we confine ourselves to the analysis of this additional ip 
dependence, and ignore other effects of J?rsoc- One ex- 
ample of the ignored effects is the -0 dependence of the 
DW width A. To be strict, A varies with -0 even when 
Hrsoc = and nonzero i?Rsoc modifies the ip de- 
pendence of A. This effect is discussed in a recent exper- 
iment^. For i?Rsoc = 0, it is commonly estimated that 
the tp dependence of A does not affect the DW motion 
significantly for small H— and/or J. We expect that at 
least for small i?RsoCj this effect is still not important. 
Below we examine the small J?rsoc regime. 



1. Bloch DW 

For a Bloch DW, the magnetization at the center of 
the DW points along the z axis and we set -0 = for this 
direction. Then Hrsoc introduces an additional Zeeman 
energy i?Rsoc = — 2A/silAxJafl cost/' to the system. 
Here, the dimensionless constant = {2TrmX/h'^)aR 
measures the strength of the RSOC. Then the total DW 
energy becomes 



Eiq, 0) = y(g, 0) - 2Msnq{H - fixJ) 
+ 2AfsfJAx(0 - aRC0Sip)J. 



(19) 



To calculate Eb in the presence of Si?, we need to calcu- 
late the shifts of the saddle point and ground state con- 
figurations due to H and J as we did in Sec. Ill Al Since 
Srsoc is independent of q, it only affects the shifts of "05 
and 0G- From 5E ^ Q, tp value of saddle(ground) point 
■05(0) for finite a/f should satisfy 



"05(0) = 0SO(GO) 



Msn 
1? 



AxJ(l sin 05(G)). (20) 



Since 05(g) - 0so(GO) < 1, sin05(G) = sin[05o(GO) + 
k^s{G) ~ 0so(GO))] may be Taylor expanded. After some 
calculation, one then finds that up to 0(5^), Eb is given 
by 



Ef 



iVa - 26qoMsniH - /3x J) (21) 
+ 2[(50o - aii{coa-ipso ~ cos'ipGo)]Ms^{XxJ) 



1 



2aR 



sin 050 sin -0, 



GO 



cos 050 cos -0GO 



(AxJ)^ 



Note that the nonadiabatic STT contribution to Eb is 
not modified by the RSOC. The RSOC effect modifies the 
adiabatic STT contribution to Eb- Since the adiabatic 
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STT contribution is dependent on the 4' dependence of 
the disorder potential, the RSOC effect is also dependent 
on the -0 dependence of the disorder potential. When the 
'0-dependence of the disorder potential energy is absent, 



St/^o = 0, one finds 

Eb ^Vo - 25qoMsn{H - PxJ) 



(22) 



Note that the result does not depend on a/? . When vZ^ ^ 
but (JV'o = (also = "050 = 0), one finds 

Eb «Vo - 26qoMsniH - px-J) (23) 
- {Msn)^i^Z^{XxJ)^ (1 - anMsniyfXxJ) , 



where v,^ — v^^ 



Note that the leading effect of 



the RSOC is to introduce a correction term proportional 
to (kjiiXxJ^ ■ Oil the other hand, when vZ^ = but 
(5-00 7^ 0, one finds 

Eb »Fo - 2SqoMsn{H - PxJ) (24) 
+ 2[(5'0o - a_R(cosV'so - cos0Go)]Msf^(AxJ) 
+ MI{1^ujZ,'^{H - pxJf 

- aR{MsVtfv^'^{\xJ? 2 (sin^so - sinV-co) 

-Ms^Vq'^XxJ (COS-0SO - cos-0Go)] • 

Again the RSOC modifies the adiabatic STT effect. Note 
that all terms containing ur are proportional to either 
sin -0GO — sin ipso or cos 0go — cos V'so , both of which van- 
ish upon the average over many potential wells. 

2. Ned DW 

For a Neel DW, the magnetization at the center of the 
DW points along the x axis and we set = for this 
direction. Then the Zeeman energy i^Rsoc due to i?Rsoc 
becomes i^Rsoc — — 2Msr2AxJa_R sin?/;. Following the 
same procedure as above, one obtains the energy barrier 
up to 0{aR), 



Ep 



iVo - 2Msmq{H - fixJ) (25) 
+ 2[5%j}Q - (5/?(sin0so - sin0Go)]A/sf^(Ax-'") 



1 



25fl 



cos V'SO cos 0GO 



'{Msn\xJ)ai 



sin V'so sin 0go 



Similarly to the Bloch DW, the RSOC effect on the 
Neel DW appears through the adiabatic STT contribu- 
tion to Eb , and is dependent on the ip dependence of the 
disorder potential energy. When the -i/^-dependcncc of the 



disorder potential energy is absent, = Sipo = 0, one 
finds that Eq. ([^5]) becomes equivalent to Eq. (P^ . Note 
again, that the result does not depend on ur. When 
ly^^ ^ but (500 = (also 0go = 0so = 0), one finds 



Eb ^Vo - 2Msn6q{H - f3xJ) 

-{MsnfyZ^{l-2aR){\xJf. 



(26) 



(27) 



Note that the leading effect of the RSOC is to introduce a 
correction term proportional to aR{\xJY ■ On the other 
hand, when vZ^ = but Sijja ^ 0, one finds 

Eb ^Vo - 2Msndq{H - ^xJ) (28) 
+ 2[(50o - ai?,(sin0so - sin0Go)]Msf^(Ax>/) 

+ aRiMs^fv^'^iXxJf [2(cosV.so - cos0go) 
+MsVtvQ^\xJ (sin0so - sin^co)] ■ 

Again the RSOC modifies the adiabatic STT effect. Note 
that all terms containing olr in Eq. (|28p vanish upon the 
averaging over many potential wells. 



III. DW CREEP IN 2D 

When the thickness or the width of a magnetic 
nanowire is larger than the collective length Lcoi , the sys- 
tem is not a ID problem any more. Here we assume that 
the width is sufficiently larger than Lcoi and the thick- 
ness is sufficiently smaller than Lcoi, so that the system 
becomes a 2D problem. In the 2D regime, the DW con- 
figuration can be described by two functions, q{z) and 
■0(z), where z denotes the coordinates along the nanowire 
width direction. In this Section, we examine the DW 
creep in this 2D regime. We find that the -^-dependence 
of the disorder potential energy again plays important 
roles, similarly to the ID case. Previous studies^^i^ of 
the DW creep motion have ignored the 0-dependence of 
the disorder potential energy. 



A. Effective energy barrier 

When the nanowire width w is larger than Lcoi, an 
entire DW line does not move simultaneously. Instead, 
a DW motion consists of a segment-by-scgmcnt motion 
of DW segments of finite lengths. In this situation, the 
thermally activated DW motion involves DW segments 
of all possible segment lengths and the DW creep veloc- 
ity is governed by the bottleneck process with the largest 
energy barrieri^^. Hence, the effective energy barrier 
Eg°°^{H, J) for the DW creep motion, which determines 
the DW velocity v{H, J) oc exp[-£;^''°°P(i?, J)/kBT], be- 
comes the maximum value of Eb{L) with respect to L, 
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Fig 3: (Color online) (a) Schematic plot of the coordinates 
system. (b) Schematic illustration of a DW segment of 
length L, which makes a thermally-assisted transition from 
the original local minimum configuration [{gm(2)},{'i/'m(-z)}] 
to another local minimum configuration [{qm,2{z)},{tpm,2{z)}] 
through the saddle configuration [{gs(2:)},{i/'s(2:)}]. The up- 
per (lower) panel shows the change of {q{z)} {{tpiz)}) during 
the transition. The areas of the grey regions in upper and 
lower panels correspond to Uq{L)XL and u^{L)\L, respec- 
tively. 



where Eb{L) represents the effective energy barrier for a 
DW segment of length L. 

Figure E^a) depicts schematically the DW configura- 
tion in the 2D system. According to Ref. H^, the eS'ec- 
tive energy i?[{(7(z)}, {?/'(2:)}] of a given DW configura- 
tion [{q{z)}, {ij{z)}] is given by 



E 




dq 



^ COS 2il} + Vdu 



dz 



(29) 



-Mstf{H - PxJ)q + Mstf^XxJj , 

where J measures the DW elasticity and K± denotes the 
DW anisotropji^i. In Eq. (P^ . the first, second, and 
third terms represent the DW elastic energy, the DW 
anisotropy energy, and the disorder potential energy, re- 
spectively. The last term in Eq. denotes the efi^ec- 
tive energy due to the adiabatic STT and the second 
last term denotes the combined effect of the Zeeman en- 
ergy due to H and the effective energy due to the non- 
adiabatic STT. One remark is in order. As in the case 
of one-dimensional DW depinning in Sec. Ill Al the ef- 
fective energy E in Eq. is a multi-valued function 
since E[{q{z)},{^{z)}] ^ E[{qiz)},{^{z) + 27r}] while 
two configurations [{g(z)}, {?/»(z)}] and [{q{z)},{tp{z) + 
2tt}] are physically identical. Nevertheless this multi- 
valuedness problem does not cause any ambiguity in the 



determination of Eb{L) in Eq. ((36|) since tp is strictly 
confined to values much smaller than 7r/4 in our study. 

In general, Vdis will depend on both q and -0, Vdis = 
Vdis{q{z),'ilj{z), z). Later we find that ijj dependence 
can generate interesting contributions, just as it did 
in the ID system. For definiteness of the illustration, 
we consider a particular type of the ip dependence of 
Vdis, arising from the position- by-position fiuctuation of 
K±. Then the fiuctuating part 6K± generates the con- 
tribution — {SK±/ Ah) 008211; to Vdis- This fluctuation 
can arise, for instance, from position-by-position fluc- 
tuations of the saturation magnetization and nanowire 
cross-section. For simplicity of the analysis, we ignore 
the fluctuating part SK± for a while and consider it in 
the later part of the analysis. 

To evaluate E, it is useful to decompose it into 
two pieces E[{q{z)},{yj{z)}] = E.Mz)}] + E4{iP{z)}], 
where 



E„ 



E,i 



dz 

T 



dz 

T 



j_ 

2h \ Xdz 



Vdis 



(30) 



-Mstf{H - l3xJ)q 



2h V d 



El. 

Ah 



COS 2ip 



(31) 



+MstfijjXxJ 



As outlined above, to evaluate Eg°'^^, we first need to 
calculate the effective energy barrier Eb{L) that a DW 
segment of finite length L experiencesi^. Suppose a DW 
segment of length L {Q < z < L) makes a thermally- 
assisted transition from one local minimum configura- 
tion [{qm{z)}T{'4'm{z)}] of the effective energy E to an- 
other local minimum configuration [{qm^2{z)}, {i^m,2{z)Y\ 
through the saddle point configuration [{'Zs(^)}, {'0s('Z)}] 
[Fig. EJb)]. These three configurations differ in the 
range Q < z < L but are essentially the same in 
the range z < and z > L since only the DW 
segment of length L makes a thermally-assisted tran- 
sition. Then the energy barrier becomes Eb{L) = 
E[{qs{z)}, {%l;s{z)}] - E[{qm{z)}, {ipmiz)}] and it can be 
decomposed into two pieces, £^<j[{(7s(z)}] — Eq[{qm{z)}] 
and £'^[{'0s(2;)}] - E^[{'ip,n{z)}]. 



1. q degree of freedom 



First, we evaluate Eq[{qs{z)}] - Eq[{qm{z)}]. The 



last term in Eq. (|30l) gives rise to the contribution 
-MstfiH - l3xJ)uq{L)L, where u,(L) = J^"- f[q,{z) - 
Qm{z)]/L[Fig.^h)] measures the typical value of the dif- 
ference qs{z) — qm{z) in the region < z < L. Since 
qsiz) — qm{z) ~ for z < and z > L, it is evident that 
Uq{L) is a growing function of L (Fig. [3]). According to 
the theory of interfaces in disordered mediaAi where the 
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disorder and the elastic energy compete, Uq{L) grows as 



a power law Ug{L) 



{L/Lc)^ where u^q is a char- 



acteristic scaling constant, C is the wandering exponent 
and Lc is the Larkin lengtb^ii^iii. For DWs formed in 
metallic ferromagnetic films, ( = 2/3^iiiii^iiiiS. To find 
out the total contribution of all three terms in Eq. ((30|) 
to Eq[{qsiz)}] - Eg{{qm{z)\\, we note that Eg{{q{z)\\ 
has the same form as the DW free energy for the purely 
field-driven DW motion. This problem has been ana- 
lyzed in Ref. [13, and we borrow the calculation result 
of Ref. [13 to obtain the characteristic L dependence of 
i?,[te(z)}]-i?,[{g,„(z)}]. 



(32) 



MstfiH - pxJ)ug{L)L, 



where the DW energy density e^i = J/2hX^. Here the 
first term includes the combined contribution of the first 
two terms in Eq. pOp . 



2. ^/j degree of freedom 

Next, we evaluate E.^,[{ips{z)}] — E.^[{il},ri{z)}]. For a 
purely field-driven DW motion, ip degree of freedom does 
not play any role for the DW creep motion if the system 
is in the regime below the Walker breakdown (the same 
holds for the DW depinning in ID systems as well, see 
Sec. ini. Then, E^,{->ps{z)}] — E^[{ijjrniz)}] is essentially 
2gj-Q8jiiii6ji7i52 ^jj^^g central task is to determine the 
effect of J on this difference. An injection of J induces an 
excitation of ip. Since the DW anisotropy {—K±cos2ip) 
favors ip = 0, the growth of "0 is strongly suppressed 
when K± is large, which is the conventional situations in 
metallic ferromagnetic systems (in ferromagnetic semi- 
conductors, K± is usually much smaller and this may not 
be the case). Then we can fairly assume that jV'l < 7r/4 
during the DW motion. This assumption is valid even 
when spatial fluctuations of K±^ exist, provided that the 
magnitude of the K± fluctuations is sufficiently smaller 
than the spatial average of K± . Under this assumption, 
cos 2ip in Eq. ([3T|) may be Taylor expanded to obtain 



E,i, — — 



J_ f drp\ K±iq,z) 2 
2h\dz) 2h 



(33) 



^Mstf^XxJ 



dz K^{q, z) 



A 



4ft 



where the position dependence of K± is made mani- 
fest. The last term of Eq. ((33)) can be absorbed to 
VdisiQ, z) in Eq to define a new effective disorder poten- 
tial V^^^iq, z), FJ|™((Z, z) = VdUq, z) - K^{q, z)/4ft. As 
long as K^{q,z) has the same statistical properties as 
Vdis('7, 2), the L-dependence oi Eq\{qs{z)}] - Eq[{qm{z)}] 
in Eq. (|32p remains essentially the same. Then we may 
forget about the last term of E^ in Eq. and consider 
only the first three terms. 



To obtain the L-dependence of £'i^[{'0s(z)}] — 
E^[{ipm(,zy\], we first examine characteristics of the sad- 
dle and minimum configurations. At these configura- 
tions, SE^j/S'tp = 0. Thus ^ps and Tpm satisfy 



ft dz^ 



TTTT + '^ ^ + MstfXx-J - 0, (34) 



where q in K±{q, z) denotes qm{z) and qs{z), respectively, 
for tp = ipm{z) and ips{z)- We analyze Eq. (p4)) under the 
boundary condition, tpmiz) — tpsiz) « for 2 < and 
z > L. Equation ([34| is solved firstly for J = 0. Note 
that Eq. has the same structure as the Schrodinger 
equationi^ -{h^ /2m)d^^ /dz^ + [U{z) - E]^ ^ for 
a quantum mechanical particle of the mass m subject 
to the potential energy U{z) with the total energy E. 
In this analogy, K ^{q{z) , z) / h corresponds to the differ- 
ence U{z) — E. In quantum mechanics, it is well-known 
that when the total energy E is smaller than the po- 
tential energy U{z), the solution ^(z) is a sum of two 
exponentially growing functions; one growing as z be- 
comes more positive and the other growing as z becomes 
more negative. For both exponentially growing functions, 
the rate of the exponential growth is roughly given by 
y/2m[U{z) — E]/h . This knowledge of the Schrodinger 
equation is directly applicable to Eq. (p4| since K± stays 
positive for all z. This analogy implies that small change 
in Kj_ within < z < L causes an exponentially large 
change in tp at the boundaries z = and L (large L limit 
is important for the DW creep motion). Combined with 
the boundary condition, and recalling that Eq. ()34|) is a 
linear homogeneous equation, we then find that both tps 
and 1pm should be essentially zero. All other solutions 
of Eq. (p4|) cannot satisfy the boundary condition and 
moreover violate the assumption \tp\ <C 7r/4 due to their 
exponential growth. 

Next, one considers nonzero J. Since Eq. (p4)) is then 
a linear inhomogeneous differential equation, its gen- 
eral solution is a sum of the general homogeneous so- 
lution for J = and a particular solution for J ^ 0. 
Due to the exponential growth, the general homogeneous 
solution should be set to zero again and we need to 
find one nonzero particular solution, which is consistent 
with the boundary condition and satisfies the assump- 
tion \tjj\ ^ 7t/4. While the exact form of the particular 
solution is difficult to obtain, it is evident from the struc- 
ture of the linear differential equation Eq. ([M)) that the 
particular solution ip should be proportional to J. Thus, 
ips (X J and 1pm oc J. As for the L-dependence of tps and 
Tpm, it is evident that they cannot grow as a power law of 
L since they arc strictly bounded below 7r/4. Thus we ob- 
tain ips oc L'^J and 'ipm oc L'^J. The proportionality fac- 
tors of ips and tpm are different since K±{q, z) in Eq. ((M]) 
amoimts to K^{qs{z), z) and K^{qm{z), z), and they are 
different. Then it is straightforward to verify that in the 
evaluation of E^[{ips{z)}]—E.^,[{'4)m{z)}], each of the first 
three terms in Eq. (j33p generates the contribution pro- 
portional to LJ^ for ip = Ips and ip = ipm- Then the char- 
acteristic L dependence of E"^ [{V's(^)}] — E^[{fpm{z)}] 
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may be expressed as 
E4{yjs{z)}] - E^[{^^{z)}] ^ MstfXxJu^{L)L, (35) 

where u^,{L) ~ J^^ ^[ipsiz) — i^miz)]/ L scales as with 
the proportionality constant scaling as J^. We remark 
that for E^[{^s{z)}] — £'^[{7/'m(z)}] to have a nonzero 
value, it is crucial to take into account the g-dependent 
fluctuation of K±. Without it, ips = 4'm and u^{L) = 
since K±{qs, z) = K±{qm, z) and both ips and ipm satisfy 
the exactly same equation [Eq. (p4| ]. 

One remark is in order. In Ref. [2^, E^[{il)s{z)Y\ — 
Ejjj[{4'm{zY\] was evaluated to be proportional to JL, 
which is different from our evaluation result J^L. This 
difference stems from the fact that the thermally- 
activated transition process considered in Ref. [H is qual- 
itatively different from the transition process considered 
in our work; While |V'| is assumed to remain smaller than 
7r/4 for the transition process considered in our work, it 
is assumed in Ref. [23l that il> jumps by ~ tt for each tran- 
sition process. Such transition with the jump of ip by 

TT may be relevant for a DW motion in ferromagnetic 
semiconductors where the magnetic anisotropy is much 
smaller. 



B. Creep velocity 

The DW velocity v{H, J) in the creep regime is given 
by w oc exp(— iJ^^^'^^/ZcBT), where for given H and 

J is the maximum value of Eb{L) with respect to L. 
By combining Eqs. ([5^ and ([55)) . we obtain the effective 
energy barrier Eb{L) for the DW segment of length L. 
Its L, J and H dependence can be summarized as 

Eb{L) ^e,i{uq{L)YL-^ (36) 
- Mstf{H - l3xJ)u,{L)L 
+ MstfXxJu4L)L, 

where Uq{L) ~ Uqo{L/Lc)'' and u^,{L) — u.^qIP J . Sub- 
stituting these relations into Eq. ([55)1 leads to 

Eb[L) ^e,i^^L^^-^~Mstf(H~PxJ)^L^^' (37) 
+ MstfXxJ'^Uii^oL. 



and Eg^^^ is written as 



For metallic ferromagnets^'^ii^'^'^ with ( 
Eq. ([57)1 becomes 

Eb{L) = AL^/^ - BL^/^ + CL 



2/3, 
(38) 



where A = e^ju'^.nLr.'^^^ , B Mstf{H - PxJ)uqaL, 



and C = MgtfXxu^oJ^- The maximum energy barrier 
Eg'^"'^ is then determined by Eg'^"'^ = Eb{Lcoi), where 
the collective length Lcoi satisfies dEs I dL]^^^^ = 0. 
From Eq. psp . the collective lengtbii Lcoi is given by 



L 



col 



-3C + V9C2 + 20AB j 



-3/2 



2A 



(39) 
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3/2. 



(~2C + V9C2 + 2QAB) 
(-3C + V9C2 + 20AB)3/2 ■ 



(40) 



1. Effective magnetic field 

The effective magnetic field H* {H, J) for the DW creep 
motion is defined by the relation v{H, J) = v{H* , 0) with 
the constraint H*{H,0) = H. The effective magnetic 
field H* provides a convenient way to express the re- 
sult for v{H, J); Recalling that the DW velocity for the 
purely field-driven DW motion is giveriF^ by v{H, 0) = 
Vq exp{—KH^'^/kBT), the DW velocity for general H and 
J can be expressed as 



v{H, J) — Vq exp 



k[H*{H, J)] 
keT 



(41) 



where k is a constant independent of H and J. Thus 
the evaluation of H* {H, J) amounts to the evaluation of 
v{H, J). H*{H, J) also determines contour lines of equal 
DW velocity in the {H, J) plane. 

Since v{H, J) is determined by Eg'^'^iH, J), H*{H, J) 
can be calculated from i;^'''=P(if, J) = E''j^''''P{H* ,0). We 

define D = MstfUgoL'^^^ and e = /3x- Then Eq. ^ 
can be expressed as 



E 



where 



creep 



-(2A)3/2(20Ai:')-i/4[F(i7, J) 



-1/4 



(42) 



FiH, J) = 



-377 J2/10 + ^(377j2/10)2-h (iJ-eJ) 



-?/J2/5 + v/(377J2/10)2 + (iJ-eJ) 



(43) 

where -q = u^oXLcxi^^^st f / ^eiu^qoY''^ ■ It can be easily 
verified that F{H, J = 0) = H. Since the constants A 
and D arc independent of H and J, F{H, J) itself is the 
effective field, H* = F{H,J). One also finds that k in 
Eq. gl]) is given by k = i2/5){2A)^/^{20AD)-^/^ . In the 
limit i7, J — > 0, we expand F{H, J) to obtain 



H*{H, J) = H-eJ--nJ^^/H-eJ+-{r]JY+0 (J*^ 





(44) 



Again, as the DW depinning in ID systems (Sec. Ill C\ . 
the non-adiabatic STT (e J) acts in the exactly same way 
as the magnetic field (H). The adiabatic STT contri- 
bution [rjj'^), however, introduces the nonlinearity to 
H* and thus plays a qualitatively different role from 
the magnetic field for the creep motion. If an exper- 
iment is performed for sufficiently small H and J, so 
that the nonlinear contributions in Eq. (|44p arc negligi- 
ble, the creep motion will follow a simple scaling behav- 
ior, v{H, J) = vo exp[-K{H - ejy /kBT] with = 1/4. 
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However if H and J are not sufficiently small, the nonlin- 
ear contributions in Eq. (|44p introduce deviations from 
the simple scaling behavior and should be taken into ac- 
count in an experimental analysis. 

C. Rashba spin-orbit coupling effects 

The RSOC is generated when the inversion symmetry 
is broken^. When a current flows in a nanowirc with the 
broken inversion symmetry, the magnetization feels as if 
there is an additional magnetic field -ffRSOC, of which 
magnitude is proportional to J—. We consider the case 
where the inversion symmetry along the ij axis is broken 
and the current flows along the x direction (parallel to 
the DW motion direction). Then -ffRsoc is along the z 
direction. When the RSOC is strong, it may modify the 
nature of the DW motion qualitatively. But when the 
RSOC is weak, its effect may be accounted for perturba- 
tively. Below we assume the RSOC to be weak. Then 
its effect can be calculated in a way similar to the ID 
case discussed in Sec. Ill Dl Again the RSOC effect varies 
depending on the magnetic anisotropy and the DW struc- 
ture. We confine ourselves to nanowires with the PMA 
and consider two types of DW structure; Bloch DW and 
Neel DW. 



The second integral of Eq. (|47p can be absorbed to Vdis 
in Eq in Eq. (|30|) . and we may concentrate on the first 
integral of Eq. (|47|) . Note that the contribution from 
the RSOC (oc Jauip'^) has the same structure as the 
DW anisotropy contribution (oc K±ip'^). Thus the main 
effect of the RSOC is to renormalize K± to £,K±, where 
^ = 1 -I- 2hMstfXxJaR/K±. Since ^ — 1 oc ArJ, it is 
safe to assume ^ — 1 <C 1 in the creep regime where J is 
small. 

For ctR = 0, it has been demonstrated [Eq. (|35p ] that 
E.^[{'tps{z)}] ~ Ejplltpjriiz)}] = MstfXxJu^iL)L, where 
= u^oJ- For nonzero an, the RSOC effect will appear 
through the rcnormalization of u^q. It is reasonable to 
expect that the rcnornialized u^o depends on ^ — 1 in a 
nonsingular way. Then we may Taylor expand u^q with 
^— 1 as a small variable, and express the renormalized u^q 
as UTi,o[l+^jiJ+0{J'^)]. Although the exact evaluation of 
7/{ is difficult, it is evident that it should be proportional 

to OtR. 

In the presence of the RSOC, the energy barrier Eb{L) 
in Eq. (f57|) is modified to 

Eb{L) =,,,^l2C-i - Mstf{H M^L^^' (48) 
+ MstfXxJ^u^o{l+lRj)L. 



1. Bloch DW 

The magnetization direction at the center of the Bloch 
DW points along the z direction. In the convention where 
■0 = for the Bloch DW, an additional Zeeman energy 
^Rsoc due to the RSOC effect becomes 

f dz 

Ersoc "2 / —MstfXxJaRCOs^, (45) 

where cur is the dimensionless RSOC coefficient aR = 
{2TTmX/h?)aR. Depending on the sign of ctR, the RSOC 
may enhance or suppress possible deviations from -0 = 0. 
Since i^Rsoc depends only on ip, it may be included in 



Ej,. Then Eq. (|3Tj) is modified to 



E,b = — 



dz 

X 



J_ /90 

2h V 



— — cos zip 
Ah 



(46) 



+MstfXxJiip - 2aR cos-0) 
For \^\ <C 7r/4, it reduces to 



, dz 
E.„=l- 



J_ fd^pV K^{q,z) , 
2h\dz) 2h 



(47) 



+MstfXxJ{i^ + aRi:'^ 



dz 

T 



K^{q,z) 

Ah 



2MstfXxJaR 



Since Eq. (|48|) has the same structure as Eq. (|37|) except 
that the last term of Eq. pS)) is multiplied by the extra 
factor (1 + 7_R J), the energy barrier E'^^'^'^^ for the creep 
motion can be obtained straightforwardly from Eq. <\A2\ . 
For metallic ferromagnets with C — 2/3, the effective field 
for the Bloch DW in the presence of the RSOC is given 
by the equation identical to Eq. except ry is now 
replaced by 77(1 -I- 7^; J). The leading correction due to 
the RSOC(oc 7_r) appears in terms of cubic and higher 
orders of J and thus, we conclude that the RSOC does 
not modify the DW creep motion qualitatively in small 
J regime. 



2. Neel DW 

The magnetization direction at the center of the Neel 
DW points along the nanowirc direction (x direction). 
In the convention where ■0 = for this direction, an 
additional Zeeman energy i?Rsoc due to the RSOC effect 
becomes 



RSOC 



dz 

'x 



■AIstfXxJctR sin-i/j. 



(49) 



Note that this equation differs from Eq. (|45p (sin0 vs. 
cos ip) since -0 = represents the different directions (z 
vs. x) in the two cases. For |0| ^ 7r/4, E^ in Eq. (|3ip is 
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modified to 



dz 

T 



2h\dz 2h 



(50) 



dz K±{q, z) 



A 



Note that an appears only in the second line, which ac- 
counts for the adiabatic STT effect. It is then evident 
that the RSOC renormalizes the adiabatic STT effect by 
the renormalization factor (1 — 2a b). 

With this knowledge, the energy barrier E'^^^'^'^ can be 
obtained in a straightforward way. For metallic ferro- 
magnets with C = 2/3, the effective field for the Neel 
DW is given by the equation identical to Eq. ([^^ except 
replacing -q by 77(1 + 2aii). Note that the correction by 
nonzero an again appears in rather high order terms in 
J. Thus we conclude that the RSOC does not modify 
the creep motion of the Neel DW qualitatively. 



IV. CONCLUSION 

Magnetic DW motion in a nanowire was examined in 
the weak driving force regime with particular attention to 
metallic ferromagnets, where the DW anisotropy is very 
large. Effects of the magnetic field, the adiabatic STT, 
and the nonadiabatic STT on the DW motion were ex- 
amined under the assumption that the amplitude of the 
tilting angle dynamics is much smaller than 27r. To be 
more specific, we examined two phenomena, the DW de- 
pinning from a single potential well in ID systems, and 
the DW creep motion through a disordered potential pro- 
file in 2D systems. 

The analysis on the ID depinning becomes relevant 
when both the width and the thickness of a nanowire are 
smaller than the collective length Lcoi- The nonadiabatic 
STT has the same effect as the magnetic field, and to- 
gether, they generate the leading order contribution to 
the depinning rate. We found that the way that the adi- 
abatic STT affects the DW depinning depends on the 
nature of disorders. In particular, it was demonstrated 
that in certain types of disorders, the conventional ways 
to determine the nonadiabaticity parameter /3 result in 



incorrect values. Possible ways to avoid the incorrect 
evaluation have been proposed. 

The analysis on the 2D creep motion becomes relevant 
when the width of a nanowire is larger than Lco\ while 
the thickness remains smaller than Lcoi- A thermally- 
assisted DW velocity is determined by the energy bar- 
rier i?^'"'''' between two spatially adjacent local minimum 
configurations in the DW energy profile. The contribu- 
tion of the non-adiabatic STT (oc /3J) to E'^^'^'^^ is the 
same as that of the magnetic field. The role of the adia- 
batic STT, however, is qualitatively different from those 
of the non-adiabatic STT and the magnetic field. Effi- 
ciencies of driving forces (magnetic field and current) are 
described in terms of the total effective magnetic field. 
Both the magnetic field and the non-adiabatic STT gen- 
erate linear contributions to the total effective magnetic 
field, implying that the purely field-driven and purely 
current-driven DW creep motions belong to the same uni- 
versality class. The adiabatic STT, on the other hand, 
generates J-quadratic or higher order contributions to 
the total effective magnetic field, and thus its contribu- 
tions constitute the next leading order contributions. Al- 
though these contributions are irrelevant in the vanishing 
driving force limit, their effects may need to be taken into 
account in practical scaling analysis since experiments 
are always carried out at small but finite driving force 
strength. 

Effects of the Rashba spin-orbit coupling (RSOC) on 
the DW depinning in ID systems and on the DW creep 
in 2D systems are also discussed. For a Bloch wall in 
a nanowire with the PMA, the RSOC effect appears in 
terms of cubic and higher orders of J in the effective 
energy barrier. For a Neel wall in a nanowire with the 
PMA, the RSOC affects the effective energy barrier in a 
way similar to the adiabatic STT. Thus, its contribution 
to the energy barrier appears in quadratic and higher 
orders of J. 
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